Superlens in the time domain 
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It has been predicted theoretically and demonstrated experimentally that a planar slab supporting 
surface plasmons or surface phonon polaritons can behave as a super lens. However, the resolution 
is limited by the losses of the slab. In this letter, we point out that the resolution limit imposed by 
losses can be overcome by using time-dependent illumination. 
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An important property of surface waves is their ability 
to focus the energy at subwavelength scales. This prop- 
erty is a consequence of the existence of surface waves 
with wavevectors K parallel to the interface larger than 
2tt /A. For an air-metal interface, the well-known disper- 
sion relation K = (w/c)\/e/(e + 1) shows that when e 
approaches —1, the wavevector modulus diverges. The 
availability of arbitrary large wavevectors is the basis of 
the plasmonic super lens proposed by Pendry (TJ [5] and 
demonstrated by several groups [3H5] ■ While experiments 
with subwavelength resolution have been reported, the 
improvement appears to be limited by the presence of 
losses as originally predicted |T] and analysed in several 
contributions [6 -8 J . Further developments have been fo- 
cussed on other schemes such as the hyperlens [HI ITU] . 
In this letter, we show that using a time-dependent in- 
cident pulse allows improving the spatial resolution. We 
report a theoretical and numerical analysis of the imaging 
properties of a SiC planar slab using time-dependent illu- 
mination. We find that the resolution can be improved if 
the pulse duration is shorter or on the order of magnitude 
than the surface wave decay time. 

At first glance, the idea of improving spatial resolution 
by using a time-dependent illuminating field seems du- 
bious as a time-dependent field is a linear superposition 
of monochromatic fields, each of them being subject to 
a resolution limit due to losses. Let us first present two 
qualitative arguments suggesting that spatial resolution 
might indeed depend on the time shape of an incident 
pulse. It is well-known that the resolution limit can be 
analysed in terms of the exponential decay exp(— Kd) of 
the evanescent waves of a field at a distance d from the 
object [TTJ [12]. In his original proposal of a perfect lens, 
Pendry pointed out that in the electrostatic regime (for 
any large wavector K), a silver slab with thickness h has 
a p-polarized transmission connecting the amplitudes at 
both interfaces given by 



t(K,u) = 



4e exp(-Kh) 
(e + l) 2 - (e - l) 2 cxp(-2Kh) ' 



(1) 



The key point is that for e(oS) + 1 = 0, the transmis- 
sion factor takes the value exp(Kh) so that the expo- 
nential decay exp(— Kh) appearing in the numerator is 



replaced by an enhancement factor exp(Kh) that com- 
pensates the decay of evanescent waves in the vacuum. It 
is useful to introduce the transmission factor t(K, oj) = 
t(K,ui) exp(-Kh) which takes value 1 for e(u>) + 1 = 0. 
As discussed in the original paper and confirmed by latter 
analysis and experiments, the presence of losses prevents 
to obtain e(w) = — 1 so that there is a practical limit to 
the resolution. 

There is another point of view that is often used to 
explain the resolution limit due to losses. When plotting 
the dispersion relation of a surface plasmon propagating 
along a metal-vacuum interface, assuming that the fre- 
quency is real, the relation K = (uj/c)y/e/(e + 1) yields a 
complex wavevector. When plotting the dispersion rela- 
tion in the plane (Re(K), u)), a backbending is observed 
(see Fig.l). Hence, the dispersion relation shows a max- 
imum value of the wavector K max corresponding to the 
turning point. From Fourier analysis, a resolution limit 
Ax = ir/K max is expected. 

The point that we want to make here is that these ar- 
guments are implicitly assuming that the measurement 
is performed in stationary regime with a harmonic field. 
Let us revisit these arguments allowing w to be a com- 
plex frequency. Equation e(w) + 1 = can now be solved 
exactly if we seek a complex solution w sw . Of course, a 
wave varying as exp(— iu> sw t) with Im(oj aw ) < diverges 
for negative times so that it is not an acceptable solution 
for all times. However, we may consider approaching this 
function in a limited time window. Improving the resolu- 
tion will depend on our ability to approach this ideal so- 
lution. Obviously, such a solution has a time-dependent 
amplitude which indicates that a time- dependent illumi- 
nation has to be used. Finally, this first remark also pro- 
vides the relevant time scale for the time-dependent illu- 
mination: the surface wave decay time l/Im(u> sw ). 

Let us now revisit the discussion of the dispersion re- 
lation for lossy materials. While solving the equation 
K = {uj/c)y/e/{e + 1) with a real frequency and a com- 
plex K yields a dispersion relation with a backbending, 
it turns out that by solving the same equation with a 
complex frequency and a real wavevector K yields a dis- 
persion relation without backbending as shown in Fig.l. 
Hence, the dispersion relation for complex frequencies 
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displays an asymptote with arbitrary large wavevectors. 
There is a price to pay. The surface waves with very large 
wavevectors which may contribute to a highly localized 
field are very short-lived modes. Note also that they have 
an almost zero group velocity. This is not a drawback for 
the superlens imaging process which is essentially an elec- 
trostatic image. For plasmons, the typical time scale is on 
the order of 10 fs whereas for surface phonon polaritons, 
decay times are on the order of 1 ps. Again, we see that 
it is possible to recover large wavevectors and therefore 
improve the resolution by working with time-dependent 
fields. Let us note in passing that the dispersion rela- 
tion with an asymptote not only leads to the prediction 
of large spatial resolution but also to large electromag- 
netic density of states close to the interface. This large 
density of states can be independently calculated in the 
direct space and has a divergence |13j . Two well-known 
physical consequences of this divergence are the extremly 
short decay time of a two-level system close to an inter- 
face [U] and the divergence of the thermal radiation den- 
sity of energy close to the interface [TS]. The existence of 
two dispersion relations which can both being meaning- 
ful depending on the context was recognized in the early 
days of surface plasmons [16] . A detailed discussion can 
be found in ref. 1171 . 
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FIG. 1: Dispersion relation of a surface phonon polariton at 
the interface SiC/air displaying the real part of uj versus the 
real part of K. Green/dashed : the wavector is chosen to be 
complex, Blue/plain : the circular frequency is chosen to be 
complex. 

From the previous discussion, we are led to revisit the 
resolution of a super lens consisting of a thin symmetric 
slab supporting surface waves. A number of conditions 
need to be fulfilled in order to achieve super-resolution: 
i) The carrier frequency has to be chosen close to the 
surface wave frequency tu sw satisfying e(u sw ) + 1 = in 
order to ensure an efficient excitation of the surface waves 
with large wavevector ii) the contribution of the transient 
surface wave needs to be the leading contribution to the 



field. If a stationary field builds up, it may overwhelm the 
transient contribution. Since the transient surface waves 
decay exponentially with a typical time scale given by 
1/Im{uj sw ), we can separate them from the incident field 
using a shorter pulse so that the resonantly excited sur- 
face waves are the only contributor to the near field after 
the end of the incident pulse. This amounts to perform 
a time filtering of the large K which contribute to the 
large resolution. This is the major difference with previ- 
ous work that considered a time analysis of the transient 
regime [HUB] or a relatively long pulse illumination[TU]. 
For plasmons, the second condition implies pulses shorter 
than 10 fs. For surface phonon polaritons, pulses dura- 
tion need to be on the order of a picosecond. 

We now show a numerical simulation of the image of 
a time-dependent point-like dipole p/(i)<5(r — r'). Since 
we are working in the non-retarded approximation, the 
field spatial structure can be derived using electrostat- 
ics. We derive the potential (f)(r,r ,t) created at r by 
a time-dependent charge qoS(t — t') located at r . The 
potential <f) p (r, r , t) produced by a dipole p can then be 
derived using p ■ V r „ { ^ I%ro ' t - ) . The electric field is then 
given by E(r, t) = — V r </> p (r, r , t). The system stud- 
ied here is a SiC slab with thickness d in a vacuum as 
depicted in Fig. 2. The dielectric constant is given by 
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cm -1 , ojlo = 969 cm' 1 , and T = 4.76 cm -1 . For 
this material, the condition e(uj) + 1 = is satisfied 
for lu sw with Re(u> sw ) — 947.7cm -1 . The correspond- 
ing decay time is \/2Im{uj sw ) = l.lps. The potential 
produced by a harmonic unit charge at r = (x, y, z) with 

P 



\/ x 2 + y 2 in a vacuum is given by: 



(r,w) 



r dK J (K P )cM-K[z 
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47T60 jp* + [ Z - (z + 2h)} 2 ' 



(z + 2h)]) 
(2) 



In the presence of a slab, the potential produced on 
the other side of the slab is found by introducing the 
frequency-dependent transmission factor: 



(r, W ) =SgMfZ°dKi(K )U )J (Kp) 
exp(-K[z- Oo + 2h)}). 



(3) 



The value of this filter is plotted in Fig. 3 for three 
different frequencies as a function of the wavector K. For 
the complex frequency lu sw , the filter does not depend on 
K . For real frequencies, the filter is band limited. Out of 
resonance, the filter decays very fast. As expected, the 
cut-off frequency is larger for Re(oj sw ) — 947.7cm -1 . 

In practice, we use a time-dependent signal composed 
of real frequencies. Hence the key issue is the design 
of the best possible time shape of the pulse. Here, we 
have not attempted to optimize the pulse shape. For 
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FIG. 2: Scheme of the planar lens with thickness h. The 
image is located at z — zq + 2h. 
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FIG. 4: Incident field H(-t) exp(-iw, 
system mode with complex frequency. 
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times t > 0. In order to compare the resolution at dif- 
ferent times t > 0, we have compensated the exponential 
decay of the amplitude. It is clearly seen that the width 
at half-maximum of the intensity distribution decreases 
as time increases showing that the resolution is increased. 
It is seen that for a 16 ps delay, the width is reduced to 
60% of its value for a monochromatic illumination. Of 
course, the price to pay for this resolution enhancement 
is the exponential decay of the amplitude. Note however 
that the incident field is extinguished so that detecting 
the signal is only limited by the signal to noise ratio. 



FIG. 3: Transmission factor as a function of the wavevector 
for three different frequencies. An ideal superlens is charac- 
terized by a transmission factor t(K,u) = 1 at the optimum 
frequency. In the presence of losses, we obtain the best be- 
haviour at w = Re(u} sw ). 

the sake of illustration, in what follows, we consider a 
stationary monochromatic field which is turned off at t — 
0. The surface wave field then naturally decays with a 
frequency close to u) sw and we study the field for t > 
0. The carrier frequency is taken to be the real part of 
the surface wave resonant frequency ui' sw — Re[uj sw ] for 
a single interface so that the incident field is given by 

We are now ready to evaluate the field in time domain. 
The potential is given by: 

/oo 
— 0(r,w)exp(-*cjt) (4) 

The results are shown in Fig. 5. We consider a dipolc 
oriented perpendicularly to the interface. We plot the 
square of the z-component of the electric field averaged 
over a few cycles. The intensity is shown at different 
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FIG. 5: Normalized profile of the intensity at the image loca- 
tion. A 60% reduction of the width of the pulse is observed 
at t = 16ps. The amplitudes habe been normalized in order 
to compare the widths. 

So far, we have introduced the idea that time- 
dependent illumination allows exciting surface waves 
with no limitation regarding the wavevector. Yet, we 
have not analysed the detailed mechanism. Analysing 
the behaviour of the planar lens requires to examine the 
image formation in terms of spatial frequencies filter- 
ing in the time domain. To proceed, we assume that 
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the time dependent source qof(t) has been induced by 
an incident electric field E(i) = E /(£). We introduce 
the frequency spectrum of the time dependent charge 
q(t) = q J f(uj) exp(— iujt)duj /2ir . The time-dependent 
potential transmitted by the superlens is thus given by : 
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dKJ (Kp) exp(-K[z - {z + 2ft)]) 



x / t(K,Lj)f(uj)exp(—iu}t) 
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(5) 



The last integral appears to be a time- dependent filter 
Il(K,t). Hence, a proper choice of f(u>) may allow tai- 
loring the spatial frequency filter. The filter H(K, t) can 
be evaluated analytically for a given pulse. The trans- 
mission factor t{K, uS) of a symmetric slab has poles as- 
sociated with the symmetric and antisymmetric surface 
modes. Their excitation generates two contributions os- 
cillating at the surface waves frequencies. A beating is 
observed as already discussed in refs [51 [TB] . Here, we see 
that the choice of the carrier frequency and the enveloppe 
of the pulse provides an additional degree of freedom to 
control the modes excitation. The structure of the filter 
is shown in Fig. 6 for the step function illumination. It 
is seen that the filter has a broader bandwidth than the 
stationary filter at lo sw . The closed-form expression can 
be found in ref.|20 . 
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FIG. 6: Spatial frequency filter H(K, t) at different times. It 
is seen that the bandwidth increases. This filter depends on 
the choice of the incident field. 

In summary, we have put forward the idea that spa- 
tial resolution using lossy surface waves can be improved 
by using time-dependent illumination. The fundamen- 
tal idea illustrated here is that when dealing with time- 
dependent fields, it is relevant to think in term of modes 
with real wavevector and complex frequency (i.e. lim- 



ited lifetime). The associated dispersion relation does 
not present any spatial frequency cut-off. A direct sim- 
ulation confirms these ideas and shows a 60% reduction 
of the width of the image of a point source. Further op- 
timization of the time shape of the incident pulse should 
allow improving further the resolution. Yet, the general 
trend is that the increase of the spatial bandwidth comes 
at the cost of an exponential decay of the amplitude of 
the field. Hence, although the resolution is not theoreti- 
cally limited with this time-dependent scheme, the signal 
to noise ratio introduces a severe limitation. 
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